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^ , 

^ We study the electronic structure and correlations in the geometrically frustrated two dimen- 

, sional checkerboard lattice. In the large U limit considered here we start from an extended Hubbard 

04 ' model of spinless fermions at half-filling. We investigate the model within two distinct Green's func- 

tion approaches: In the first approach a single-site representation decoupling scheme is used that 
^ , includes the effect of nearest neighbor charge fluctuations. In the second approach a cluster repre- 

' sentation leading to a 'multiorbital' model is investigated which includes intra-cluster correlations 

' exactly and those between clusters on a mean field basis. It is demonstrated that with increasing 

nearest-neighbor Coulomb interaction V both approaches lead to a metal-insulator transition with 
^SJ ' an associated 'Mott-Hubbard' like gap caused by V. Within the single site approach we also explore 

the possibility of charge order. Furthermore we investigate the evolution of the quasiparticle bands 
as funtion of V. 
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I. INTRODUCTION 



i 

i' ' To treat the problem of a metal-insulator transition (MIT) driven by electron-electron interactions in lattice models 
'■^ , with a fractional electron site occupancy, both the on-site and long-range Coulomb repulsions are equally important 
Ch ' Q- A lattice with geometrical frustration imposes additional complications and may produce new features of MIT. 
O ] The well known example is the Verwey transition Isl in magnetite {Fe^iOij at Ty ~ 120K. Magnetite has a spinel 
structure and A sites are occupied regularly with Fe^'^ ions. On the B sites forming a pyrochlore lattice, the Fe 
ions are in the mixed valent state of Fe^+ : Fe^+ = 1:1. Therefore, in the high-temperature charge disordered state, 
T > Ty, a band structure calculation would predict a quarter-filled conduction band, i.e., a metallic state with one 
^ ' itinerant electron per two sites. Below the Curie temperature Tc ~ SSOiiT of magnetite the itinerant electrons are 
(almost) fully ferromagnetically polarized. Therefore they can be treated as a system of interacting spinless fermions. 
In terms of spinless fermions, the charge occupancy of one electron per two sites corresponds to a half band-filling. 
, Below the Verwey temperature Ty, the charge degrees of freedom arc ordered and an insulating state occurs. As 
^-H ' was first pointed out by Anderson , properties of the Verwey transition can not be understood without taking into 
' account the geometrical frustration and the resulting huge degeneracy of the ground-state charge configurations in 
the pyrochlore lattice. Two more recent examples of electronic structure changes connected to charge ordering in 
spinels are found in AIV2O4 with temperature decreasing at ambient pressure 0,"^ and in LiV20i under the external 
pressure {t^ . An interesting observation is that charge ordering in these three spinel systems is accompanied with a 
lattice structural change, i.e., the system tries to avoid the geometrical frustration of the pyrochlore structure. At the 
same time, no MIT and charge ordering were observed in LiV204 at ambient pressure down to very low temperatures, 
but instead the metallic compound LiV20i exhibits a heavy-fermion behavior below 30K 
Q Most of the work until now has been devoted to an understanding of the magnetic properties of geometrically 

O frustrated lattices such as pyrochlore structure, since in the presence of antiferromagnetic interaction frustration acts 
^ ' against a conventional long-range order and may stabilize a spin- liquid state 0. Charge degrees of freedom have 
• ^ in contrast been studied much less lllll. In fact, charge ordering in geometrically frustrated systems has been an 
intriguing and unsettled problems |12|. For example, numerical diagonalizations of a Hamiltonian for spinless fermion 
^ \ with strong nearest-neighbor repulsion on a checkerboard lattice have given evidence that at half filling (nurnber of 
■ - - ' spinless fermion equals one half the number of sites) the ground state is two- fold degenerate but a liquid d3] . This 
must be kept in mind when that system is treated within different approximations leading to a charge ordered ground 
state. In other words, the observed charge order in frustrated structures could crucially depend on associated lattice 
deformations, i.e., the involvement of lattice degree of freedom. 

The extended Hubbard model on non-frustrated lattices has been extensively investigated in one dimension at 
quarter- or half-filling 15], for two-leg ladders at quarter-filling 16j, for two-dimensional square lattices at half- 
filling^ | and in the limit of infinite dimensions at quarter- |18J or half- filling jl^. A variety of techniques, such 
as Hartree-Fock approximation, pertubation theory, dynamical mean-field theory, as well as numerical methods, e.g., 
quantum Monte Carlo and density-matrix renormalization group have been employed. However all these investigations 
were based on the non-frustrated lattice. 
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In the present paper, we apply a Green's function approach to study the possible phase transitions of a half-filled 
spinless fermion model , i.e., a system with one electron per two sites, on the frustrated checkerboard lattice (see 
the inset of Fig. 1). The lattice can be viewed as a two-dimensional projection of a pyrochlore lattice. The spinless 
fermion model arises naturally for ferromagnetic materials in which one of the spin-split bands is completely occupied 
or completely empty as in magnetite ,20 1 . Also this model can be viewed as a quarter- filled extended Hubbard model 
in the large U limit. Then double occupancy of a site is forbidden and the nearest-neighbour Coulomb repulsion V 
plays a crucial role. Because for (spinless) half filling every second site is unoccupied on the average one would naively 
expect a metallic state. Our main goal in this work is to show that inter-site correlations V can lead to a MI transition 
even for a case with less than one electron per site. We are using two different methods to study this model. 

Firstly we employ a single site approach within a Hartree-Fock as well as Hubbard-I type approximations fSect. lTTIl . 
By using the 'Hubbbard-I'-type decoupling scheme which includes the effect of nearest-neighbor charge fluctuations, 
we flnd that with increasing value of V first a metal- insulator transition occurs with a gap in the excitation spectrum, 
while at even larger values of V charge ordering appears. This is opposite to the result in the simple Hartree-Fock 
approximation, i.e. inter-site correlations favor the MI transition and suppress the CO. This observation may indicate 
that CO is indeed not present for a rigid checkerboard or pyrochlore lattice for any V/t ratio (it certainly is not in 
the limit V/t^ oo). Indeed in the compounds AIV2O4 and LiV20i (under pressure) where CO has been found it is 
accompanied by a lattice distortion. 

Secondly we start from a cluster representation of the model where the intra-cluster correlations are taken into 
account exactly and the inter-cluster terms are treated in Hartree Fock approximation (Sect. !!!!)) . This transformation 
leads to an effective multi-orbital extended Hubbard model. Again we find a M-I transition at a value similar to the 
first approach. In view of the suggestion above we do not consider the possibility of CO in this case, although within 
the Hartree Fock approximation for the inter-cluster interactions it would presumably be present in the ground state. 

Furthermore, we investigated the evolution of quasiparticle bands in the various phases in the single-site and cluster 
approaches. Following the selfconsistently determined chemical potential and the formation of interband gaps allows 
to determine the critical value for the MI transition. 



where V is the Coulomb repulsion between nearest neighbors denoted by {ij). We refer to this spinless Hamiltonian 
as the t — V model. 

First, let us briefly review previous work on the t — V model at half- filling and the closely related quarter-filled 
extended Hubbard model in the large U limit for different lattices. In one dimension, the t — V model can be mapped 
onto an anisotropic Heisenberg model[23| and solved exactly via Bethe ansatz. In this case, one flnds a gapless 
metallic phase for V < 2t and a gapped charge-ordered state for larger values of V. For t — V model on a two-leg 
ladder, it was found from renomalization group calculations that contrary to a single chain, the ladder becomes a 
Mott insulator for arbitrarily small repulsive interactions y|22j|. By using the density matrix renormalization group 
method, Vojtail6j et. al. studied the extended Hubbard model for two-leg ladders in the large U limit. They found 
that the charge-ordered phase vanishes for V < 2.5t but claimed that there will be a charge gap for all values of V/t. 
For a square lattice, McKenzie|23| et al. argued from slave-boson theory that the insulating phase with charge order 
is destroyed below a critical value V of order t and the system becomes metallic. 

Here we will study the t — V model on the checkerboard lattice at half-flUing (one electron per two sites) with the 
aim to flnd out about possible MI and CO transitions in this frustrated lattice which is a 2D model for the pyrochlore 
lattice. At flrst we will employ the Green's function method and various decoupling schemes within the single-site 
representation. There are two sites per unit cell. For convenience, we rewrite the Hamiltonian (1) in the following 
form 



II. SINGLE SITE APPROXIMATION 



The Hamiltonian for the spinless fermion model is given by [l, 




(1) 



H = Ht + Hv,MF + 



V 



(2) 



where 
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,1 + Cl-y,l + Cl+x,2) + 4,1 ("^1,2 + Cl-|-x,2 + Q-ya) + h.C. 
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FIG. 1: Charge order parameter n as a function of V/t within a mean-field approximation. Dotted fine indicates the critical 
point of the metal-insulator transition. Inset is an illustration of the checkerboard lattice with (right) and without (left) charge 
order. The 1-th unit cell is indicated by dotted lines. The wave vector of the staggered CO is Q = (0, 0) due to the two sites 
per unit cell. Here charge order occurs before the metal-insulator transition takes place. 



Hv,MF = V ^[ni,2 ((ni_x-y,i + ni_y,i -|- ni+x,2)) + (^1,2) (ni_x-y,i + ni_y,i + ni+x,2) 
1 

((ni,2 + ni-Kx,2 + ni-y,i)) + ("1,2 + ni-Kx,2 + m-y,i) 

- (m,2) ((ni_x-y,i + rai-y,i + m-Fx,2)) - (ni,i) ((ni,2 + ni+x,2 + ni_y,i))], (4) 

Hy = ^^^[^1,2 ((5ni_x-y,i + fci-y,i + fei-i-x,2) + ^m,i (^ni,2 + ^ni-i-x,2 + (5ni_y,i)]. (5) 
1 

Here we have introduced a charge fluctuation operator Srii^i = ni^i — {ni^i) on sites i = 1, 2 of the 1-th cefl. is the 
kinetic energy term, Hv,mf is mean-field part of the interaction term while H'y is the residual interaction part. In 
the following we only consider the simplest charge ordered pattern, namely a staggered checkerboard pattern with 
wave vector Q= (0,0) illustrated in the inset of Fig. 1. The average charge density on diff'erent sites in the cell is 
('^i.i) = \ — where the order parameter n means a charge disproportionation within each unit cell. 

The electron propagation is described by a retarded Green's function (for simplicity the conventional symbol R is 
omitted) 

G,,,(l-l',a;) = ((c.,|4,))^= j ((cH(t)|4.))e-*rfi. (6) 

where 

' ;cH(i)l4,)) = -ie (i) {{cu {t) , 4,- (0)}) . (7) 

and {, } denotes the anticommutator. The above (2 x 2)-matrix Green's function must satisfy the equation 

uj{{A\B))^ = {{A,B}) + {{[A,H]\B))^. (8) 
By introducing the Fourier transformation 

(k, ^) = ^ E (1 - 1', c) , (9) 

k 

the equation (8) can be now written explicitly as 

(wT-hA(k)) G(k,w) =T-^yr(k,w). (lo) 



with 



f -{3-2n)V + 2tcosky 4ie-'^^ cos % cos ^ \ 
A(k)= \k^+k' ,1, , (11) 

\ 4te'^^ cos 1^ cos ^ -{3 + 2n) V + 2t cos J 
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FIG. 2: Quasiparticle dispersion within mean-field theory for different value of V/t. The chemical potential /i is fixed at zero 
energy, (a) is the noninteracting case without charge order. These two bands touch at (tt, tt) and the upper band is fiat, (b) 
is in the metallic phase with small charge order. The two bands still touch while the upper band become dispersive due to the 
inequivalence of diagonal hopping, (c) is exactly at the critical point. The two bands touch at (0, tt). (d) is in the insulating 
phase. The two bands seperate. 

where 1 is the unit matrix and F (k, uj) is the Fourier transformation of the (2x2 )-matrix F^.j (1 — 1', uj) {i, j = 1,2) 
which is defined as 

F,,,(l-l',c^) = //cu ^ 5n(„.„.)|4A\ . (12) 

\ \ (n.n.)e(l,0 / / ^ 

Here the summation is over six nearest-neighbor sites (n.n.) surrounding the site (1, i) 

A. Mean-field approximation 

In the first step we may decouple 

F,.,(1-1',C^)« / '5"(n.n.))((cMl4,,))^=0, (13) 

\(n.n.)6(l,0 / 

which leads to a mean-field approximation and the quasiparticle dispersions are given by 

E^/t = -{cosk^ + cos ky - 3V/t) 

±[4(1 -I- cosfca;)(l -I- cosfcj^) -I- {2nV/t - cosfc^, -I- coskyf]'^/'^ . (14) 

Now the intra-cell charge disproportionation n has to be determined from the following self-consistent equation 

1 {2nV/t — cos kx + cos ky) 

" ~ 2iV ^ [4(1 + cosA::,)(l + cos fey) + {2nV/t - cos/c^ + cos ky)^]^/^ ' ^ ^ 

The calculated phase diagram is shown in Fig. 1. In the present mean- field approximation (F = 0) the metallic phase 
is already charge ordered before a charge-transfer-type metal-insulator phase transition takes place at V = 2.8t. In 
the non-interacting case V = (Fig. (2a)), the upper band is flat while the lower band is dispersive. It touches the 
flat band at (tt, tt). As shown in Fig. (2b, 2c) with increasing intersite repulsion V the flat band becomes increasingly 
dispersive with increasing charge order and the touching point moves from (tt, tt) towards (0, tt) . The dispersion of 
the previously flat band is due to the inequivalence of diagonal hopping (see inset of Fig. 1) induced by charge order. 
Finally, for V > 2.8t the two bands separate as shown in Fig. (2d). In the half-flUed case, the lower band is fully 
occupied and the upper band is empty resulting in a charge-transfer-type insulator. 
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B. 'Hubbard I'- approximation for the inter-site correlations 

The neglect of correlations overestimates the tendency to CO symmetry breaking. Therefore, in the strongly 
correlated case {V ^ t), the mean- field results are unreliable and it is necessary to consider equations of motion of 
higher-order Green's functions T.^ j (1 — 1', w) which can be written as 

(^-V {n(n.n.))\ ((cl^ ^ (^"(n.n.) l4 j ) ) 

(n.n.)e(l,j) / ^ ^ (n.n.)e(l,i) ' ' u> 

(n.n.)e(l,i) (n.n.)G(l,i) 

Y (("-!:«) "("»•"■>) ((c(«.n.)l4,, 
(n.ri.)e(l,i) 

X! X! (l - '^(n.n.)',(n.n.)) (('^("•"■)"^"'("-"-)l4j, 

(n.n.)e(l,i) (n.n.)'G(l,i) 

H H ((ci.»cfn.«.)^(»-"-)'l4j] 

(n.n.)^(l,'i) (n.n.)'e(n.n.) 

+^ H H (l-'^(n.«.)',(M)) ((ci,<cJ,.„.)'^(»-"-)l4,j^ 

(n.n.)6(l,i) (n.n.)'£(ji.n.) 

+ f E -^"(n.n.) I l4,, )) ' (16) 

\ \ \in.n.)e(U) J 11^ 

Here we adopt the following approximations, which consist in an extension of the on-site Hubbard I decoupling 
scheme p3| to spinless fermion on a checkerboard lattice with intersite Coulomb interaction: 

In equation the fourth term involves both the nearest- and next-nearest-neighbor charge correlations and we 

neglect the latter. When treating the last term in equation (|16ll we also neglect the more distant charge correlations 
and approximate as 

( (cu E '^"("■"•) l4,w ) 

\ \ V(n.n.)e(l,l) / 11^ 

)2 
) ((^''114,,))^ 
/ 

|~6n2 + 2/<5n,,2 ^ '^^(«.".)) ] ((cm|4,,))^ • (18) 

\ (ri.n.)e(L2) / / 

Finally the equation of motion for the higher-order Green's function Tij (k, uj) can be written as 

{uA + A (k)) r (k, cj) = S (k) + M (k) G (k, oj) . (19) 



where 

^fcx+fcji J,^ ^^^^ 

8ni ^^"^ cos ^ cos ^ VKi 



M(k) = ( _ , -^"'^"^ ^.^^ cos -cos- I ^ (20) 
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FIG. 3: Results within Hubbard I-like approximation, (a) charge order (CO) parameter n versus V/t. The dotted Une indicates 
the phase boundary between the metalhc and insulating phase without charge ordering. Contrast to the mean-field approxima- 
tion, the charge-ordered phase transition is of first order and occurs after a Mott-Hubbard-like metal-insulator transition takes 
place. The dashed line shows the quasiparticle gap versus V/t in both insulating phases, (b) Hopping amplitude (solid line) 



from site 2 to site 1 



C1V12 



and nearest-neighbor charge fluctuation correlations Snu ^ 



(n.ri.)e(l,i) 



V/t. In the insulating phase without charge order, the expectation value (^c[^c\2 
keeps constant and maximum. Note that Snu „ )g(i i) ^^(n.n.) 
between ordered and disordered state versus V/t. For V > Vc2, the charge-ordered state has lower energy. 



Sn(^n,n.) (dashed line) versus 

equal to while Snu J2{n.n.)e{hi) ^"(n.u.) 
is site-independent for all value of V/t. (c) Energy difference 




2 (^cuqLyi^ cos ky 

-j. \ ^ -\- k'l 



4 ( cue 2 ) e 



4(cii4 



cos ^ cos ^ 



COS -Y cos -J- 



2 (C12CJ+X2) cos h 



(21) 



K, 



2 ( Jnu 



5n 



{n.n. 



(n.n.)e(l,i) 



(22) 



By solving equation (|19() with respect to F and after substituting F into equation (|10|l one obtains the final solution for 
the Green's function G. For the electron concentration 1/2 (one electron per unit cell), the chemical potential fi, the 
charge disproportionation n and the hopping amplitute ^cnCjj^ are calculated from the following set of self-consistent 
equations: 



1 



—y 

N/2 ^ 



lm[Gii(k,u;)+G22 (k,c^)]. 



(23) 
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n = tttt: / dijJ f | Im 



Gil (k,w) - G22 (k,Lj) 



(24) 



Ciic|2 



N/2 ^ 



duj 



ImGi2(k, w). 



(25) 



The higher-order Green's function Ti j (k, oj) can be easily derived from equation pU|) . Therefore the nearest- neighbor 
charge fluctuation correlation (5n\^i ^ ^^.j (5n(„ „ ) 



can be calculated as 




Sn 



(n.n.) 



N/2 



duj 



IvaVi^i (k, w) . 



(26) 



Furthermore the hopping amplitude \ c\ic\_^^j and yCi2C][_^x2/ '^^'^ determined from Gn (k, cj) and G22 (k, a;) 

respectively. 

The one-particle retarded Green's function exhibits four poles given by the roots of D = where D is defined as 



D=(uo-Et 



V^KiK2 + IGn^V^t^ (1 + cosfc^) (1 + cosfc^) 



-V^Ki[{lu ~ {3~2n)V + 2t cos ky) - (1 + cosfca;) (1 + cosfcj^)] 
-V'^K2[iuj - {3 + 2n)V + 2t cos k^f - (1 + cos k^) (1 + cos ky)] 



(27) 



As a consequence, four bands are obtained. For = 0, two branches of momentum-dependent dispersions are obtained 
from equation (|27|l : E^\v=o — —2t{coskx + cos ky + 1) and E^\y=o = 2t as expected. When t = 0, two poles ui — 2V 
and uj = AV are obtained from above equations, which describes one hole or one particle excitation. 

In Figure (3a) we show the phase diagram within the previously used decoupling scheme. First, a metal-insulator 
transition occurs at Vd = 4.86t and then charge ordering appears at Vc2 = 6.47i. In the interval V > Vd and 

V < Vc2 a gap opens and increases with increasing value of V. At Vc2, it drops to a smaller value. For V > Vc2, the 
gap increases again monotonously. In the charge ordered phase, there exist always two self-consistent solutions, i.e., 
a charge ordered and a disordered one. By comparing the energy of these two states, (see Fig. 3 (c)), we find for 

V > Vc2 that the charge ordered state has lower energy. The phase transition is of first order. Fig. 3 (b) shows that 

at the metal-insulator transition the hopping amplitude ^ciic|2^ vanishes. 

In the following we want to explain this feature in more detail. For that purpose we study the poles of the retarded 
Green's function as a function of increasing value of V/t. Consider first the trivial case of = (see Fig. 2 (a)). 
The states in the dispersive band are of bonding and in the flat band of antibonding character. As seen in Fig. 
4 (a), (b) and (c) these bands split into four when < V < Vd and the Hubbard-I like approximation is made. 
The retarded Green's function has therefore four poles for each k vector. In the regime < < Vd, states in the 

antibonding flat band become more and more occupied as V increases. Therefore the expectation value (c|, C12 



decreases since an equal occupational probability of a bonding and antibonding state implies that ^Cj-^ci2y = 0. This 

is the case in the regime Vd < V < Vc2 shown in Figs. 4 (d), (e). Note that a correlation gap has opened above 
the two occupied lower bands which implies a Mott- Hubbard type insulating state and the nearest-neighbor charge 
fluctuation correlations bnn 'Yliin n )e(i i) ^^{n.n.) remain constant. When V > Vc2, then due to symmetry breaking one 

can no longer distinguish between bonding and antibonding state and (^cl^cu'^ 7^ (see Figs. 4 (f)). As V continues to 
increase, charge order become more and more pronounced and (c^-^cu/ decreases again monotonously. This explains 



the behavior of (cl-^cuj shown in Fig. 3. The inflexion point at V/t — 2.2 shown in Fig. 3 (b) is due to the lower 
flat band crossing the inflexion pionts (maximum of DOS) of the upper dispersive band shown in Fig. 4 (b). For 



Sn 



{n.n. 



which 



V > Vc2, CO strongly suppresses the nearest-neighbor charge fluctuation correlations Snu J2(n n )g(i %) ' 
is site- independent for all value of V/t as expected. 

In this section we have studied the half-filled (one electron per two sites) spinless t — V model on a checkerboard 
lattice within mean field theory and a 'Hubbard-I type' approach. In the former approximation which overestimates 
the tendency to symmetry breaking, CO appears before the system becomes an insulator. By using a Green's function 
approach and a decoupling scheme that includes the effect of nearest-neighbor charge fluctuations, however, it was 
shown that fir.st a MI transition into an insulating state without CO takes place and only for larger V eventually 
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FIG. 4: evolution of the quasiparticle spectrum within a Hubbard I-like approximation for different value of V/t. The chemical 
potential fi is fixed at zero energy. V/t = is already shown in Fig. 2 (a). Figs, (a), (b) and (c) are in the metallic phase 
(0 < < Vci = 4.86t). Figs, (d) and (e) are in the insulating phase without CO symmetry breaking (Vci < V < Vc2 = 6A7t). 
Reminiscent of the upper and lower Hubbard band splitting, the metal-insulator phase transition at Vci is of Mott-Hubbard 
type. Fig. (f) is in the insulating phase with CO {V > Vc^)- 



CO appears. It is not entirely clear, however, how realistic the appearance of CO is. In the hmit (V/t — + oo) when 
the hopping vanishes, the ground state for half filling is macroscopically degenerate 0, 0| and no CO is present 
[23. It is therefore possible that CO obtained here for the rigid checkerboard lattice above Vc2 niay be due to the 
employed Hubbard-I type approximation. In fact, the charge ordered state obtained here is one of the macroscopically 
degenerate states which must obey the tetrahedron rule (Anderson rule) . If there would exist corresponding lattice 
distortion, for example, compressed along the diagonal direction the system would select such a CO state shown 
in the inset of Fig. 1 out of the macroscopically degenerate states. Indeed in some 3D pyrochlore compounds which 
exhibit CO as mentioned in the introduction, it is always accompanied by symmetry lowering lattice distortions which 
remove the frustration by introducing inequivalent bond lengths in the tetrahedrons of the corner sharing lattice. 



III. CLUSTER APPROXIMATION 



We want to supplement the previous calculation based on a Hubbard I type of approximation by another one where 
the strong electron correlations are treated exactly within a cluster but within mean-field approximation outside the 
cluster. Within this scheme we want to determine the critical interaction Yd at which a gap opens in the excitation 
spectrum when the case of half-filling is considered. We do not care here about a possible charge order at large value 
of V because, as pointed out before, that may turn out to be an artefact of the involved approximations. 

We divide the checkerboard lattice into sublattice A and B of the plaquette so that each sublattice contains 
units where N is the number of sites. Fig. 5 shows black (basic) clusters linked accross white square. The Hamiltonian 
(1) , H = Ht + Hy, is decomposed into 

Ht = ijf "*™) + ; Hv = + (28) 

Here each of the terms, fj^"^*^'^^ and jg ^ sum over I = 1, ...,N/A decoupled basic clusters while Tj^f*"*'^'") and 



FIG. 5: Fragment of the checkerboard lattice. Only four neighboring basic clusters connected by elementary translation vectors 
n, T2 and ra are shown. Within each basic cluster the lattice sites are denoted by 1,...,4. 



^(inter)^^^ due to the intercluster coupling. Specifically the decoupled clusters are treated by the Hamiltonian 



^i^a) (J) ^ ^ /_^^t^^^.^ ^ y^nan^ . (30) 

Here i, j = 1, ...,4 denote the sites within a cluster with lattice vector 1 and (0) is a (4 x 4)-matrix with elements 
^ij (0) = Ij if * 7^ j-i ^^^d Mii (0) 0. Furthermore nn = c\^Ci\. Thus, fj^"]^"-") (i) represents an intracluster 
Hamiltonian to be solved below. 

To write down in a systematic way the intercluster coupling given by iff "'^'■^ and iJ^;"*^'"^ , we refer to Fig. 5 where 
a fragment of the checkerboard lattice is shown with the origin located at the center of a basic cluster. There are 
eight neighboring basic clusters connected to a given one with the translation vectors {t } — Ti, T2,...,T8. Three of 
them are depicted in Fig. 5. The intercluster electron hopping and Coulomb repulsion terms can be written as 

^i.nter) ^ ^ ^ ^ (^) ct ,^^c,, (31) 

1 {T} »J 



1 {T} 



(32) 



where the eight (4 x 4)-matrices My (r ) are specified by referring to Fig. 5. Consider first the matrix Mij (xi) 
that connects two neighboring clusters by Ti-translation. According to Fig. 5, only two individual two-site bonds 
contribute to this connection (reading from the right to the left): (zj) = (14), (23). Therefore, we define Mi4(xi) = 
M23 (ti) — 1 and (ti) = otherwise. In a similar way, one finds M13 (T2) — 1, Afi2 (T3) = M43 (ra) = 1 and 
Mij (T2) = Mij (ra) = otherwise. With this procedure, the other five matrices Mij (r ) can be easily found as well. 

Now we solve the intra-cluster eigenvalue problem by diagonalizing the Hamiltonian }jj^"]^^°-'> (]^ from (|30() separately 
in each n-particle (n=0,l,...,4) sector 



(33) 



The empty state \ >i with energy -E^^i = defines the cluster vacuum state |0 >i. There are four singly occupied 
cluster states >i with v = 1, ...,4 which are found to be 



*W >,=V/3,4|0>i^/t,|0>i 



(34) 



where P^i is a (4 x 4)-matrix 



^111 1 \ 

1-11 -1 

V2 -V2 

V V2 -V2 J 



(35) 
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and the corresponding eigenvalues are 

^ = -3t; Ei^'> = ) = ^ = t (36) 

The eigenvectors (|34|l are basis vectors of different irreducible representations of C^h point group: the case = 1 
describes the lowest energy fully symmetric (og) solution, 1^ — 2 belongs to bg and i' — 3, A to representation. 
Below we will refer to these 1-particle solutions as cluster 'orbitals'. According to (|34|) . for a given cluster 1 the 
transformation from the original site operators cj^j (i=l,...,4) to the cluster 'orbital' operators fl^ is given by the 
matrix /3^i defined in (|35|l . Note that the /^j (/^i) operators anticommute. The cluster states of higher occupancy 
71=2,3,4 can be understood as a result of successive filling of cluster 'orbitals'. For instance, for n—2 one obtains six 
eigenstates (the cluster index 1 is dropped): 

\^?l3 >= fUlflfl /]/! |0 >, 1*1^6 >= fill fill M\Q > (37) 
with corresponding energies 

E^L = -2i + V, , = 2t + V (38) 



In the cluster 'orbital' representation, the terms ff^^™*'"'^) and in can be now written as follows 

^(™*.a) ^ ^ l-itfljn + t ^ fl,U,] , (39) 



1 \ v=2 J 

4'"*™^-yEE E (40) 

^From the result (|4U|I . one can see that ffj^"*'''') jg i^j^g Hubbard term in an effective, 'multi-orbital' electronic model. 
Such a model is derived below by adding to H39I) and H40() the intercluster hopping ^^(™*'^'') Coulomb jj^"-*'^^^ 

terms and by using for the latter a mean-field approximation, i.e., n^^i+xTiji ~ (j^i.i+x) 'T'ji+'T'i,i+T ('^i.i+x) 
The approximated term Hy^^p reads 

h^mP -ImvY^jz = i (^^) ^ E E -^1' (41) 

1 i=l 1 iy=l 



where the last equality is due to Y^^nn = X^i/^i'i) ^"^^ i^c) = (^i/i) is the average cluster occupancy. In the 
cluster 'orbital' representation, the hopping term takes the following transparent form 

^i.nter) ^ ^ ^ ^ T^^, ) /t^^^/^„. (42) 

1 {T} vy 

The (4 X 4)-matrices T,^^/ (r ) are related to the Mij (r) by a rotation 

T,,, {T)^~tY, (r ) , (43) 

because = f^v'j- Finally, by collecting the contributions H39|) - (|42|l . we obtain an effective 'multi-orbital' 

Hubbard-like Hamiltonian 

Ht-v = ijf -f -f + irj^^t/F ■ (44) 

The effective Hamiltonian (I44II has a lower symmetry compared to that of the original Hamiltonian (1). This may 
lead to an artificial low-symmetry ground-state solution, for instance, to a long-range charge/bond ordering. Here we 
suggest, however, that the geometrical frustration of the checkerboard lattice prevents this kind of long-range ordering, 
as was discussed already in the beginning of this section. Therefore we supplement the effective Hamiltonian (|44f) with 
additional restrictions, which prevent an artificial low-symmetry solution of (|44|l . We require: (a) a homogeneous 
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electron distribution over the lattice sites, yc|ciy = n, with n being the electron concentration, and (b) a bond- 
independence of hopping amplitudes, ^c|cj^ ~ k, where a pair of sites i,j denotes a bond ( there are six bonds within 

each plaquette). Note that fc is a real n-dependent quantity to be calculated self-consistently. On general grounds, 
one expects that > in a metallic phase and fc=0 in an insulating phase. Both requirements (a) and (b) together 
impose restrictions on cluster 'orbital' averages 



{t.i^'^'^ 



(45) 



By using that fiuiPv'i = ^i/z^' and 



5I^''»-{o;' = 2,3,4} 



one obtains for v = v' the following relations for cluster 'orbital' occupancies 

J/iVii) = " + 3ft, (^^=1); (/,V-i)-"-fc, ('^^ 2,3,4). 
If 7^ i^', the expression H45|l leads to 

(/,V-^'i)=0, {v + v'). 



(46) 



(47) 



(48) 



For an isolated cluster with an integer electron occupancy equation H48I) is obviously fulfilled because of symmetry 
arguments. 

Based on the effective model (|44|) we calculate the electronic band structure from the Fourier transformation 
Giiyi (q, w) of the retarded matrix Green's function: 



G... (1 - 1', i - t') = ((/., {t) mU. it'))) = -^e it - t') {{.U it) , fl,, it')}) . 



(49) 



To obtain the equation of motion for G^,^' (q, uj), we use the methodpq of the two-time 'irreducible' Green's function. 
A successive differentiation of (|49f) with respect to both times t and t' with the use of properly defined projection 
procedure lead to the Dyson's equation: 



wl-17(q)-S(q,w) G(q,w) = l 



(50) 



Here 1, (q) and I](q, cj) are the unit, a frequency and a self-energy (4 x 4)-matrices, respectively. To define il 
and E explicitly, it is convenient to introduce a notion of a scalar product of two fermionic operators A and B as 
{A\B^') ~ (^{A,B^}'). In this notation, the matrix elements of are 



resulting in 



n,,, (q) = (5,,, + T,,, (q) + S,,,V 



Here the sj:,^'' are given by H36|l and 



T,,, (q) = e'l-^ T,,, (r ) . 

{T} 




(51) 



(52) 



(53) 



The frequency matrix $7 provides for a mean-field description of the electronic band structure. In order to include 
effects of electron correlations, the self-energy part S must be calculated. 
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For this purpose, the set of basis operators /^q is complemented with a new set of operators F,yq (^^ = 1, 4): 

which are orthogonal to /i/q, i.e. {^-Fl/q|/i!/q' ^ — 0- Then, the self-energy S has the form of an 'irreducible' matrix 
Green's function[2^ whose elements are defined by 



The irreducible Green's function matrix ('('^i/q|^"J'q*)\ obeys the following equation (the lower indices are omitted 



for brevity): 



We truncate this equation by neglecting in H56() the next-order 'irreducible' matrix Green's function 
. t \ \ 

i F \ — i F )) ■ For the effective 'multiorbital' model (|44|l . the truncation corresponds to the first step of 



the Hubbard-I approximation and results in the following form of S: 

E.,, (q, c^) = LT - A (q)] (F,,q|Fj, \ , (57) 

L J vv' \ I 

where the frequency matrix A (q) is given by 

A..' (q) = ^7 T^. (58) 

In the Hubbard-I approximation applied to the standard Hubbard model, the corresponding self-energy is q- 
indcpendent. We shall apply a similar approximation here by dropping the intercluster hopping term ^^("'*^'') which 

appears in the equation i Fuq— [Fuq., Ht-v]- In the following, A and E are q-independent matrices and the latter is 
of the form 

S..' (^) = <5.. . -Y ■ (59) 

u^-{^ Fui \Fl,) I {f,i\fI,) 

With the definition (|54|) . one obtains explicitly 

= VSNU.i. (60) 

where 

^iV- = TV- - (iVO , N^^N,~n,,^ "-ii- (61) 

Then the diagonal elements of the frequency matrix entering into the denominator of H59|l are found to be 

' - ^ - + (^c) + (62) 



F.AFI,) ^ (<5iVn 
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FIG. 6: Energy bands of spinless fermions on the checkerboard lattice calculated in the cluster approach for different values 
of inter-site repulsion Vjt = 0, 3 < Vcjt and Vjt = 10 > Vc/t ~ 4. The bands axe unfolded along the highly symmetrical 
directions in the square Brillouin zone; energy unit is t s 1. The electron concentration is chosen as n = i (equivalent to iN^ 
=2) and the chemical potential /^ is fixed at zero energy. 



Both averages, ((5iVf)^> and iN{ {SN^'f), are approximated in a mean-field manner. For instance, in this approx- 



{SN,f)= ^ 



imation, the mean value ( {5N^) ) which describes the intracluster charge correlations, reads 

("-I'll) (1 - ("-I'll)) • 
We summarize the results for the self-energy as follow ((nyi) = {ny)): 



'i^vu (w) = 



A. 



(63) 



AM K +4^^^^^ + ^E.,(^.)K.>(l-K.)) 



V 



(1 



l)){n^2) 



(64) 



In this expression the second term is the inter-cluster Hartree-Fock correction. The third and fourth terms are intra- 
cluster Hubbard-I type correlation corrections, where the latter appears only for the present 'multi-orbital' Hubbard 

Hamiltionian. In the common Hubbard model (with just one 'orbital') only the third correction would be present. 
To perform self-consistent band-structure calculations, a chemical potential n is introduced in a standard way. For 
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FIG. 7: Mean cluster occupancy (A^c) calculated as a function of the chemical potential /i at V/t = 10. Two wide flat regions 
with d{Nc)/d^ = at (A^c) =1,2 indicate an insulating state of the system. 

a given spinless-fermion concentration n varying within the range < n < 1, the value of is determined from the 
equation 



Here the summation is over M = L/A q- vectors in the first Brillouin zone of the square lattice formed by the basic 
clusters in the checkerboard lattice. The 'orbital' occupancies should obey the relations l|47(l . 

The most representative results of the self-consistent band structure calculations are shown in Figs. 6 and 7. In Fig. 
6, band spectra for different values of coupling V are given along highly symmetrical directions in the square Brillouin 
zone. To show the insulating gap opening in the band spectrum with increasing V, the electron concentration n=l/2 
corresponding to an integer mean cluster occupancy {Nc)=2 is chosen; the chemical potential ^ is located at zero 
energy. From the upper panel, a doubly degenerate completely flat branch on the top of two highly dispersive bands 
is seen for the non-interacting case, V=0. For finite coupling V, but less than some critical value V < Vc, these four 
branches are split and the chemical potential intersects one of the dispersive branches as shown in the middle panel 
of Fig. 6. In this weakly correlated regime of the model one expects a metallic state of the system. We found that 
starting from V = Vc ~ 4t, the spectrum is clearly split into the low- and high-energy Hubbard subbands separated 
with a gap growing as V increases above Vc- For a given concentration {Nc)=2 and V > Vc, the four lowest Hubbard 
subbands are filled completely and located below the chemical potential fi as shown for y = 10 in the lower panel of 
Fig. 6. If one includes a weak next neighbor hopping t' , the former flat energy branches acquire a small dispersion, 
the picture of the insulating gap opening still retains at slightly changed value of the critical coupling Ve- 
in this paper, however, we restrict ourselves to the study of the limit, t'=0. In this limit, to avoid discussing a 
rather special case of the partially filled upper flat band, the electron concentration is chosen to vary in the range 
n < 1/2. Within this range, for strong coupling V the model displays an insulating state at the integer cluster 
occupancy {Nc)=l as well. This can be easily seen from Fig. 7, where the calculated mean cluster occupancy (Nc) 
is depicted as a function of the chemical potential fi for V=10- The charge compressibility d (Nc) /d/i is found to be 




(65) 



where (n^) is an average 'orbital' cluster occupancy 




(66) 
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zero in a wide range of varying fi at both integer occupancies {Nc)=l,2. Near the occupancy value (iVc)=1.5, the 
sharp change of (Nc) is connected to the fact that /i intersects the flat band peculiarity in the density of states. 



IV. CONCLUSION 

In this work we have studied an extended Hubbard type spinless fermion model for the frustrated checkerboard 
lattice which is a 2D analogue for pyrochlorc or spinel-type lattices. We have studied the possibility of a metal- 
insulator phase transition as function of the inter-site Coulomb interaction and the band filling. We have used 
single-site and cluster representations of the model. In the latter the intra-cluster correlations are accounted for 
exactly by transforming to a 'cluster orbital' basis. Within a Green's function approach decouplings on Hartree- 
Fock and 'Hubbard-F type level for the inter-site Coulomb term have been employed. In both approaches we find 
a metal-insulator transition for increasing V for half filling (that is one fermion per two sites) and in the cluster 
approach even for quarter filling. The MI transition is of the Mott Hubbard-type and is associated with a gap 
opening in the quasiparticle excitations. The critical interaction for the MI transition is Vci/t=4.86 in the single 
site approach and Vc/t ~4 in the cluster approach are in reasonable agreement. We also consider the possibility of 
charge ordering within the single-site approach where we find a transition to a staggered CO state with Q = (0,0) at 
Vc2/t=6.47. This state may be a result of the approximations employed since for V/t oo CO has to vanish due to 
the macroscopic degeneracy of the ground state. Stabilisation of CO may indeed require the coupling to the lattice 
to lift this degeneracy . Further progress in understanding the nature of the isulating state caused by the inter-site 
Coulomb interaction may require the use of more advanced methods like cluster dynamical mean field methods which 
has sofar not been achieved. 

The authors are grateful to Prof. P. Fulde for suggesting the subject of studies and valuable discussions and 
comments. 
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